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Abstract: We construct the skyrmion solutions appearing in the coset spaces 
SU{N)/ SO(N) for N > 2 and compute their classical mass. For = 3, the third ho- 
motopy group n3{SU{3) / S0{3)) = Z4 implies the existence of two distinct solutions: the 
skyrmion of winding number two has spherical symmetry and is found to be the light- 
est non-trivial field configuration; the skyrmion and antiskyrmion of winding number plus 
and minus one are slightly heavier and of toroidal shape. For > 4, there is only one 
skyrmion since the third homotopy group is Z2. It is found to have spherical symmetry 
and is significantly lighter than the N = 3 solutions. 



1 Introduction 



Skyrmions [1] are topological solitons appearing in field theories where the set of mappings 
from three-dimensional space to a target space Ai is split into non-equivalent homotopy 
classes, characterised by the third homotopy group tts^M.). The lowest energy configuration 
in each equivalence class cannot be transformed continuously to the vacuum and is therefore 
stable over time. If those field configurations are degenerate in energy, they represent 
different vacua of the theory. If not, only the lightest one is a true vacuum state, the others 
being called skyrmions and representing upon quantisation distinct particle states. 

In the low-energy chiral theory of QCD, cr-models describing mesons are known to con- 
tain skyrmions, which represent the baryons of the theory [2, 3] . In this case the global sym- 
metry group consists of flavour SU (2), or can be extended to flavour SU (3). The non-trivial 
homotopy group 'it3{SU{N)) = 7L allows the existence of an infinite tower of skyrmions 
characterised by an integer charge, the baryon number. Baryons as skyrmions of a scalar 
theory of mesons have been extensively studied in the literature (see [4] for a review). 
Recent results include stabilisation of the skyrmion from a spatial extra-dimension [5, 6]. 

Many modern theories trying to solve the hierarchy problem raised by the Higgs mech- 
anism of electroweak symmetry breaking enlarge the particle content and the global sym- 
metry group of the Standard Model. Some of them predict new scalar fields — both 
elementary or appearing as fermion condensates — at low or medium energy scales, and 
share therefore similar features with the Skyrme model. This is the case in particular in 
models where the Higgs boson is a pseudo-Nambu-Goldstone boson of some strongly in- 
teracting theory, as in the composite Higgs models [7, 8], or more recently in little Higgs 
models [9~15], as well as in other strongly interacting composite models arising from a 
holographic description [16-18]. All these models describe the Nambu-Goldstone sector 
where the Higgs lives in terms of a cj-model, with various symmetry groups depending on 
the exact realisation of the model. Some of these symmetry group have a trivial topology 
and can not lead to the existence of skyrmions, but other do. 

In this paper we will focus on the case of a SU{N)/ SO{N) cr-model, which has a 
non-trivial topology as characterised by its third homotopy group [19] 



7T3iSU{N)/SO{N)) = Zp, p 



4 for TV = 3, 

(1.1) 

2 for TV > 4. 



Such a symmetry is realised with = 5 in the littlest Higgs model [11], for which skyrmions 
are present [20] and can even act as dark matter [21-23]. ^ In general, the SU{N)/SO{N) 
symmetry breaking pattern can be obtained from a strongly interacting theory with Weyl 
fermions in the adjoint representation of the gauge group [25-28]. The finiteness of the 
third homotopy groups (1.1) implies that SU{N)/ SO{N) skyrmions have quite different 
physical properties from their SU{N) cousins. Spherically symmetric skyrmion solutions 
of winding number B = 2 for = 3 and of winding number B = 1 for A^ > 4 have been 



^The non-trivial topology of the SU (5) / SO{5) coset also induces the presence of other topological defects 
in the littlest Higgs model, such as cosmic strings and Z2 monopoles [24]. 
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constructed before in ref. [29, 30] and [23] respectively. We review these two cases and 
complete the analysis by constructing a solution of winding number B = ±1 for A'^ = 3, 
whose symmetry is only axial. 

This work is organised as follows: in section two the model is introduced and general 
considerations about SU{N)/ SO{N) skyrmions are presented. In the third section we 
present a general method to construct solution of unit winding number and apply it to the 
cases N = 3 and > 4. The classical masses of the corresponding skyrmions are computed. 
The fourth section contains the construction of the N = 3 skyrmions of winding number 
two. Finally, in the fifth section, the results are summarised and further issues are shortly 
discussed. 



2 SU(N)/SO(N) skyrmions 

We consider a cr-model defined by the Lagrangian density 
C = C2 + C4, 

= ^—TrduT.d''J:'' + ^ Tr S'^^.S, S^^a^S S^S^S, S^'a^'E , (2.1) 

4 32e^ L J L J 

where S(x) is a SU{N) symmetric matrix. The model has a global SU{N) symmetry 
under which T,(x) transforms in the two- indices symmetric representation as 

S ^ U^U^, U G SU{N). (2.2) 

The vacuum state is obtained by taking a constant value for the field S(x). Here we choose 
it to be the identity matrix 

(S) = 1^. (2.3) 

This vacuum expectation value spontaneously breaks the global SU{N) symmetry down 
to SO{N), since only the SO{N) subgroup of SU{N) leaves the vacuum unchanged under 
the transformation rule (2.2). As a consequence the Nambu-Goldstone field S(x) takes its 
values in the coset SU{N)/SO{N). 

The Lagrangian (2.1) is identical to the one of the Skyrme model [1] describing the 
low-energy chiral limit of QCD, in which / is identified with the pion decay constant and 
e is a parameter depending on the high-energy behaviour of the theory, found empirically 
to be e = 5 [3]. The notable difference is that in the original Skyrme model the field 
transforms in the adjoint representation of SU{Nf)L x SU{Nf)ji, where A'^^ is the number 
of light flavours, and a global diagonal SU{Nf)v symmetry is preserved after spontaneous 
symmetry breaking. 

The presence of skyrmions in the model (2.1) is due to the fact that the third homotopy 
group TTs of the target space is non-trivial. However, this condition is not sufficient to ensure 
the existence of skyrmions. If the four-derivative term £4 were omitted, the mass and size of 
the skyrmions would shrink to zero, according to Derrick's theorem [31]: under a rescaling 
of the spacetime coordinates x — >• Ax, the energy associated with the kinetic term scales as 
£"2 — >• A £'2 and can thus be made arbitrarily small. Only the interplay of C2 and £4, whose 
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energy scales as — )• {X)~^Ei, ensures a stable size for the skyrmion. In the presence of 
both terms, the energy is bounded by a minimum i?min = 2 -y/ -£'2-^4, reached when E2 and 
E4 contribute equally. 

In this paper, for simplicity, we will always consider static solutions — moving 
skyrmions can be obtained by applying a Lorentz boost — and work in dimensionless units 
X = {fe)x. With this choice, the energy density corresponding to the Lagrangian (2.1) is 



E[T,] = ^ j d^xTi '^diT.diT,^ + ^ {^.iT.^^T.^ djT.djT) - diT,djT,^ diT,djT,^ 



(2.4) 



The integrand is then independent of the parameters / and e and can be directly computed 
from the shape of the skyrmion solution. 

The elements S(x) of the coset SU{N)/ SO{N) can be parametrised by the Cartan 
embedding [19, 23, 25], starting from a matrix <I>(x) in SU{N), as 

SU{N)^ SU{N)/SO{N), ^^^^ 
$(x) 1-^ S(a;) = ^{x)^{x)'^. 

The matrix ^(x) is only defined up to right-multiplication with an SO{N) matrix. The 
advantage of this embedding is that the winding number of a field configuration ^(x) in 
SU (N) can be simply expressed in terms of an integral over space, as^ 

Bi<^) = ^e^jfc J d^'x Ti{^^di^){^^dj^){^^dk^) G Z. (2.6) 

This integral is a topological invariant, and is additive with respect to multiplication of 
two fields $i(x) and (^2{x): 

B{^i^2) = B{^i) + B{^2) + ^ Jd'x di^u (2.7) 

where = ejjfc Tr($|5j<I>i)(<I>2<9fc*^*2)- '^^^ ^^^^ term is a surface term and therefore van- 
ishes provided that the fields are in the vacuum at the spatial boundaries, $1,2(2^ 00) = 1. 

The winding number of S (x) can then be straightforwardly identified with the winding 
number of the field ^(x) used to construct it. However, this mapping is not unique: the 
multiplication of $(x) from the right with an SO{N) matrix R{x) leaves S(x) unchanged, 
but not the winding number (2.6), which is raised or lowered as B{^R) = B{^) + B{R). 
The crucial point here is that B(R) cannot take an arbitrary integer value, since the 
winding number integral (2.6) evaluated on an element of SO{N) gives an integer multiple 
of 4 for = 3 and of 2 for > 4 [19]. Therefore the winding number of an element 
in the coset SU{N) / SO{N) is only defined modulo a factor of 4 or 2, for = 3 and 
> 4 respectively. Moreover, two field configuration Si(x) and S2(x) built respectively 
from SU{N) fields <&i(x) and $2(2^) differing by one unit of winding number cannot be 



^The notation B comes from the original Skyrme model of mesons, in which the winding number is 
identified with the conserved baryon number of the theory. 
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continuously transformed into another and thus belong to different homotopy classes. This 
is an immediate consequence of the statement made above that 7r3{SU (3) / 50(3)) = Z4 and 
7r3{SU{N)/SO{N)) = Z2 for > 4. The winding number for the field = ^>(x)^>(x)^ 
can thus be expressed unambiguously as 

= mod p, (2.8) 

where p is defined as in eq. (1.1). 

With the help of the Cartan embedding (2.5), the challenge of constructing a field 
configuration in each of the homotopy classes of the coset space reduces to finding an 
appropriate configuration $(a:) in SU{N) with the desired winding number. A naive choice 
consists in taking the lightest SU{N) skyrmion solution of a given winding number and 
building the corresponding S'C/(A^)/S'0(A^)-valued field directly out of it. For B = 1 the 
lightest solution is known to be the hedgehog configuration (3.2) with spherical symmetry 
and mass Mo = 72.9 I with our conventions [3]. The Cartan embedding of this solution 
yields a skyrmion of mass 4Mo = 291.7 | whose symmetry is however only axial. This 
is a hint that the hedgehog ansatz might not yield the lightest field configuration, as we 
shall see. For B = 2, the SU{N) skyrmion has a toroidal shape, and its mass of 139.6 ^ is 
slightly lighter than twice the mass of a single skyrmion [32] . The naive Cartan embedding 
of this B = 2 configuration gives a mass of 599.0 { for the corresponding SU{N)/SO{N) 
field. We can therefore infer the following upper bounds on the skyrmion masses: 

f f 
Mb=±i < 291.7 -, Mb=2 < 599.0 -. (2.9) 

e e 

We are actually going to show in the next sections that the true skyrmion solutions have 
masses much below these naive bounds. 

Notice finally that in the original SU (N) Skyrme model a six derivative term is often 
present as an alternative to the four-derivative term £4 or in addition to it. It contains only 
two time derivatives and preserves all the symmetries of the Skyrme Lagrangian [33, 34]: 

£6 = C6TVS^S'^, (2.10) 

where B^ is the topological current defined as 

^'^ " 2^2^^'''"'/ '^^^ Tr{<^^d„^){^^dp<^){^^d„^). (2.11) 

However, the equivalent of this current with the field S(a;) replacing $(x) vanishes identi- 
cally due to the structure of the coset, so that the Lagrangian term (2.10) is automatically 
absent of the theory. For the same reason, the usual Bogomolny bound [35] for the mass 
of a SU{N) skyrmion of winding number B, Mb > 6'7r'^\B\^, does not apply since the 
component Bq of the topological current (2.11) is zero. A crucial consequence is that the 
mass of the skyrmion is not necessarily increasing with increasing winding number. Indeed, 
we shall see in the following that for A'^ = 3 the B = 2 skyrmion is lighter than the B = ±1 
ones. 
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3 The solutions of winding number B = 1 



The skyrmion configuration of unit winding number in SU (N) has been known for a long 
time. It is obtained by using a so-called hedgehog ansatz living in a SU{2) subgroup of 
SU{N). Due to the symmetries of the original Skyrme Lagrangian, the embedding of the 
SU{2) subgroup can be chosen without loss of generality to be in the upper-left 2x2 block 
of the SU{N) matrix as 



Here F{r) is a function of the radial coordinate, Xi = Xi/r are angular coordinates, and the 
cTj are the usual Pauli matrices. The boundary condition for F at spatial infinity is fixed 
to -F(oo) = in order to recover the vacuum ^{x — )■ oo) = 1, while for definiteness at the 
origin F{0) has to be an integer multiple of vr. Chosing -F(O) = vr ensures the correct unit 
winding number. 

The hedgehog ansatz is built so that a spatial SO (3) rotation of the coordinates Xi — )• 
RijXj is equivalent to a SU{2) transformation <I>o ~^ Uq^qUq, with the equivalence between 
the two given by Rij = ^ Tr (aiUoCjUQ). Since the original Skyrme Lagrangian is symmetric 
under diagonal SU{N) transformations — hence under any SU{2) subgroup — it is also 
invariant under spatial rotations. For this reason, the skyrmion built using this hedgehog 
ansatz is said to be spherically symmetric. In general, soliton configurations with the 
highest symmetry tend to have the lowest energy: this is true not only for skyrmions, but 
also for monopoles and dyons. We want therefore to follow the same general principles to 
construct the SU{N)/ SO{N) skyrmion. 

The first important point about the SU {N) skyrmion construction is that it only makes 
use of a SU{2) subgroup of the target space. This is the minimal choice, as proven in a 
theorem due to Bott [19, 36]: the winding number integral (2.6) is actually only counting 
the number of times the field is winding around any SU{2) subgroup of SU{N),'^ which 
makes it sufficient to embed all the non-trivial components of the field ^(x) into a SU{2) 
subgroup, while the other can be taken in the vacuum. We postulate here that this rule 
can also be applied to the construction of the SU{N)/ SO{N) skyrmion, so that in general 
the field ^>(x) can be written as <I>(x) = exp[i /j(x)Tj], where the Tj are generators of a 
SU{2) subgroup of SU{N). 

However, while in the original Skyrme model the relevant subgroup could be trivially 
chosen to live in the upper- left 2x2 block of the matrix <I>, this is not the case in the 
coset, where not every choice of SU{2) subgroup is equivalent: among the diagonal SU{N) 
transformations $ — t- U^U\ only the ones satisfying C/^(S) = (S)C/-^ are symmetries of the 

^This is a direct consequence of the fact that the physical space is isomorphic to S'^ ~ SU{2) upon 
identification of the spatial infinity to a single point on the sphere. 




(3.1) 



where ^q{x) is the previously mentionned hedgehog ansatz described by 



^o{x) = exp [i F{r) Xi ai] . 



(3.2) 
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model. The form of the ansatz (3.1) is nevertheless very useful to ensure a correct winding 
number by construction, so instead of considering a different ansatz for $ we keep this 
form and work in a different basis where the vacuum expectation value is not necessarily 
diagonal. The most general form of the Cartan embedding (2.5) is S(x) = ^{x)'Eq^{x)'^ , 
where = (S) denotes the vacuum state. 

The second crucial point about the construction of the SU{N) skyrmion is to impose 
spherical symmetry. We have seen that this is possible if the transformation $o ~^ Uo^qUq 
is a symmetry of the Lagrangian. This SU (2) transformation of can be promoted to a 
SU{N) transformation by considering ^> — )■ U^U^ , with 



U 



Uo 



\ 



V 



Uo£SU{2), V £ SU{N -2). 



(3.3) 



As discussed above, this transformation is a symmetry of the Lagrangian only if it satisfies 



C/tSo = Sot/' 



(3.4) 



Since Uq and Uq live respectively on the left- and right-hand side of eq. (3.4) and cannot 
be related by multiplication with a matrix Eg for arbitrary values of Uq, the equation can 
only be fulfilled if the (A^ — 2)-dimensional matrix V is proportional to Uq. From this point, 
we have to consider two separated cases: 

1. For N = 3, the one-dimensional matrix V is fixed toV = l, hence eq. (3.4) cannot be 
satisfied and the SU{2) transformation rule of $0 cannot be promoted to a symmetry 
of the Lagrangian, independently of the basis. Thus we conclude that there exist no 
spherically symmetric ansatz of unit winding number in the coset SU{3)/ SO{3). 



2. On the contrary, for N > 4, choosing 



V 



1 



N-i 



u 



Uo 



\ 



ux 



\ 



1 



(3.5) 



the transformation $ — )• U^U'^ acts on <I>o as a diagonal SU{2) transformation 
<I>o — ^ Uq^qUq and simultaneously preserves the vacuum, so that the Lagrangian 
can be made spherically symmetric. This is the most general choice of basis, up to 
global SU {N) transformations, and we expect therefore to build the lightest skyrmion 
configuration with the help of this ansatz. 



3.1 The N > 4 skyrmion 

With the choice of basis (3.5), we have just proven that the skyrmion can be made spheri- 
cally symmetric. Indeed, an important property of the Skyrme Lagrangian (2.1) is that it 
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can be expressed in terms of the currents which take in this case a block diagonal 

form, 



\ 



\ 



V 



OAr-4 J 



(3.6) 



so that each of the two 2x2 blocks contributes equally to the energy density (2.4), which 
becomes twice the functional expressed in terms of the 5C/(2)-valued field ^o{x), 



(3.7) 



Hence the mass of the SU{N)/ SO{N) skyrmion with > 4 is exactly twice the mass of 
the original SU{2) skyrmion. The solution minimising the energy -E[$o] is obtained using 
the hedgehog ansatz (3.2), for which the winding number integral (2.6) is unity: 



B($) 



dr F' sin^ F 



sin F cos F — F 



F=0 



{31 



The energy density — explicitly independent of the angular coordinates due to the spherical 
symmetry — reads 



E[F] = iir 



f 



dr 



r^ + 2 sin^ F) F'^ + {2r^ + sin^ F) 



sin^ F' 



(3.9) 



and is minimised when F is a solution of the Euler-Lagrange equation 
(r^ + 2 sin^ F) F" + 2rF' + sin 2F ( F'^ - 1 - ^^^^ 



0. 



(3.10) 



The numerical solution of this equation and the corresponding energy density are displayed 
in Figure 1. The skyrmion mass computed as the total energy of this field configuration is 
then 

/ 



MN=i B=i = 145.8 



(3.11) 



This mass is exactly twice the mass of the SU{N) skyrmion, and only half the value of the 
upper bound (2.9) found from the naive Cartan embedding of the original solution. We do 
not have a proof that it is indeed the lightest topologically non-trivial field configuration, 
but it is the only simple embedding preserving the spherical symmetry of the Lagrangian, 
and it is much lighter than other constructions based on different choices of ansatz, as 
for example the = 3 case discussed below. Note that this solution is the one used in 
[23] with N = 5, although the choice of basis (S) = I5 makes the ansatz there look more 
complicated than the simple form (3.1). 
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01 2 _1 4560 J 2 3 4 

r r [1/fe] 

Figure 1. Left: The profile functions F and x for the N = A, B = 1 and N ~ 3, B = 2 skyrmions 
respectively. Right: the corresponding radial energy distributions dE — —AttC (right) as functions 
of the radius r in dimcnsionful units of 1/fe. 




Figure 2. The numerical solution in the case = 3, S = 1 for (left) and for the difference 
G — 6 (right) showed in the plane {x,y) = {f sin 9, r cos 6). 

3.2 The N = 3 skyrmion 

The solution of unit winding number in SU{3)/ S0{3) was not computed so far in the 
literature, and requires a more involved ansatz. The reason for this is that the spherical 
symmetry cannot be preserved, as argued before. 

In the original basis (S) = I3, the ansatz (3.1) nevertheless allows to preserve an 
axial symmetry in the Lagrangian, and we expect therefore that it yields the lowest-energy 
skyrmion solution. The most general form of $0 preserving this axial symmetry can then 
be written in analogy to the construction of the deuteron solution [32] as 

$o(x) = exp [i F{r, 9) mir, 9, (^) ai] , (3.12) 

where n is a vector of unit length, given in spherical coordinates by 

(ni, ^2, 713) = (sin G(r, 9) sin (/?, cos G(r, 0), sinG(r, 9) cos 93) . (3.13) 
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Instead of a single profile function F of the radial variable r as in the hedgehog ansatz (3.2), 
we are now left with two functions F and G of two variables r and 6. The boundary 
conditions are fixed to F{0,9) = vr, F{oo,6) = 0, G(r, 0) = and G(r, vr) = vr, so that the 
field ^{x) is well defined everywhere. With this ansatz, the winding number integral (2.6) 
becomes 

OO TT 

B(<P) = ^ j dr j dO sin^ F sin G{drFdeG - deFdrG) (3.14) 



and gives the expected value of one upon integration by parts. Similarly, a solution of 
winding number B = —1 (or equivalently B = 3) is obtained by inverting the boundary 
condition for F at the origin to F{0, 9) = — vr; the energy of this solution is identical to the 
B = 1 case, and we will therefore not discuss it further. Plugging the ansatz (3.12) into 
the energy functional (2.4), one obtains 



E[F,G] 



A-K j dr j ddsine I sin^ G ( 



,'2 TT' ^2/1 2 



sin^ 2Fsin'^G 



sm 



+ sin^F r^(l -cos^Fsin^G) + 



sin^Fsin^ 2G 



sm 



1 ^ ^( 1 . sin2 F sin2 G 
+ - sin 2F sin 2G + 4 ^ 



2 y sm 

+4 sin^ F sin^ G [drFdgG - dgFdrGf 
siiP F sin^ G(l - sin^ F sin^ G) 



drFdrG 



+ 



sm 



(a.G)V^ 



doFdoG' 



(3.15) 



The corresponding Euler-Lagrange equations for F and G are too long to be displayed here, 
but can be solved numerically using relaxation methods. The solution of these equations 
are shown on Figure 2, and the corresponding energy density on Figure 3. As one can see, 
unlike the spherical energy distribution of the previous two sections, the energy density is 
in this case located along a torus. ^ The mass of the skyrmion, obtained by integrating this 
energy density, is found to be 



M^=3,s=±i = 273.5 



/ 



(3.16) 



4 The solution of winding number B = 2 for N = Z 

The general rule stating that solutions with the highest symmetries tend to have the lowest 
mass also apply to the A'^ = 3 case, where the spherically symmetric solution of winding 
number two found by Balachandran et al. [29, 30] is the lightest known skyrmion configu- 
ration. It is constructed from the SO{2>) generators Tj G {A2, A5, A7}, where the Aj denote 

*Note that our solution resembles the numerical skyrmion solution of the 0(3) cr- model of ref. [37]. 
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iV = 3 
B = ±l 



N = 3 
B = 2 



N>4 
B = 1 



Figure 3. Left: The energy density of the = 3, i? = 1 solution in a plane containing the axis of 
symmetry of the skyrmion (located along x = 0). Right: isosurfaces of energy density dE = 100 ^ 
showing the relative size of the three solutions. 



Gell-Mann matrices. The ansatz is built as follows: 



$(a;) = exp [i a{r) Xi Tj] exp 



iX{r) 



;t - {x^Ti)'' 



(4.1) 



The two exponentials commute; the first of them belongs to 5*0 (3) and is therefore irrele- 
vant in the coset space, as can be seen from the Cartan embedding (2.5), yielding 



exp 



2ix(0 ( -1 



However, it is the first exponential which fixes the winding number (2.6): 



oo 

i?($) = — dr \a' (cos a cos X ~ 1) 
TT J 



x' sin a sin x] = ^ (sin a cos X — «) 

-' TT 



(4.2) 



(4.3) 



r=0 



since the expected value i? = 2 of the winding number is achieved by choosing the boundary 
conditions a(oo) = and a(0) = vr. Nevertheless, with this choice the first exponential 
in (4.1) is ill-defined around the origin, and the introduction of the second exponential with 
boundary conditions x(oo) = and x(0) = tt is required to cancel the angular dependence 
at the boundary point r = 0. Plugging this ansatz into the energy density (2.4), one obtains 
the following functional, also independent of the angular variables, 



E[x] = Stt 



/ 



dr 



+ 2 sin^ x) x'^ + (2r^ + sin^ x) 



which is minimised when x satisfies the Euler-Lagrange equation 

In 2 \ // 2 , f ,2 sin^ X 
-r + 2 sm X I X + g^'X + sm 2x I X - 1 ^ 



0. 



(4.4) 



(4.5) 



The numerical solution for xi^) is shown on Figure 1, together with the energy distribution 
along the radial direction. The mass of the B = 2 skyrmion is hence 

/ 



Mn- 



:3,B=2 



228.7 



(4.6) 
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B = ±l 


B = 2 


N = 3 


273.5 { 


228.7 ^ 

e 


N >4 


145.8 

e 





Table 1. Classical masses of the SU{N)/SO{N) skyrmions. 



5 Summary and conclusions 

The masses of the three different SU{N)/ SO{N) skyrmion solutions are summarised in 
Table 1. For = 3, the skyrmions of winding number plus and minus one are heavier 
than the one of winding number two, so that the mass hierarchy is inverted compared 
to SU{N) models. Although this result can be surprising, it could be expected from the 
symmetry of the solutions. Notice that the two masses are nevertheless close to each other 
— only about 20% difference — whereas in the SU{N) case the Bogomolny bound [35] 
implies that the mass of the second heaviest skyrmion — which has a toroidal shape as in 
the SU{3)/ SO{3) case — is nearly twice the mass of the lightest, spherically symmetric 
skyrmion. The > 4 solution with spherical symmetry is significantly lighter than both 
A'^ = 3 skyrmions. It is actually a rather surprising fact that for any values of A'^ > 2, there 
always exists a spherically symmetric solution, and that this solution is the lightest even 
when its winding number is not unity. 

The finiteness of the third homotopy groups for A^ = 3 and A^ > 4 ensure special 
annihilation properties for the skyrmions: in the latter case, the skyrmion is its own anti- 
skyrmion, so that any two skyrmions can annihilate into a final state with trivial topology. 
For A^ = 3, the skyrmion of winding number two has the same properties as the A^ > 4 
one; on the other hand, two toroidal skyrmions can annihilate into a state containing a 
spherical symmetric skyrmion or not, depending if they have the same or opposite wind- 
ing number. Notice also that although multiple skyrmion solutions are not favoured by 
energy considerations, no important attractive force between them is expected at large 
distances [23], so that the skyrmions appearing in theories described at intermediate ener- 
gies by a SU{N)/ SO{N) fi-model can be long-lived and hence provide viable dark matter 
candidates [21-23]. The A^ = 3 case is especially interesting since the two different kinds of 
skyrmion are both stable and can therefore account for a fraction of the currently observed 
relic density each. 

Note also that the masses computed in this work are only obtained at the classical level. 
The quantisation of the model typically makes the physical masses increase depending on 
the spin of the skyrmion [3]. The stability of the skyrmion in the quantised version of the 
theory is also questionable, but recent results proved that it is ensured in strongly coupled 
models where the symmetry breaking pattern SU{N)/ SO{N) arises from a condensate of 
fermions in the adjoint representation of a SU{N) gauge group [25-28]. Notice moreover 
that the inclusion of higher-dimensional terms in the Lagrangian (2.1) would also modify 
the skyrmion masses, but not their symmetry properties nor the way they are embedded 
in the coset. Finally, realistic models may require to gauge totally or partially the global 



- 11 - 



SU{N)/ SO{N) symmetry, which can affect the skyrmion stabihty and mass [23]. In these 
cases, a detailed study of the model is required. 
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